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Abstract. The two-dimensional Hubbard model on the square lattice is studied in the presence of lattice 
distortions in the adiabatic approximation. The self energy is computed within perturbation theory up 
to second order, which provides a way for studying the quasiparticle dispersion. We compute numerically 
the second order contribution to the self-energy using a standard Fast Fourier Transform Algorithm for 
finite sizes system. The stability of the lattice distortions is investigated and a schematic phase diagram is 
drawn. The Fermi surface is analyzed for densities close to half filling, the presence of lattice distortions 
changes some spectral properties of the model and gives an anisotropic interacting Fermi surface. The 
spectral function is calculated along several lines in momentum space and the renormalized quasiparticle 
dispersion is obtained. The behavior of the density of states is shown for different values of the intrasite 
repulsion U in the different phases. 

PACS. 71.10.Fd Lattice fermion models - 63.20.kd Phonon-electron interaction 

1 Introduction an antiferromagnetic Mott-insulator. The competition be- 

tween the kinetic and Coulomb terms gives rise to strong 

rrii . 1- ■ iTTuuj 11 electron-elcctron correlations, 

ihe two dimensional Hubbard model has been usually as- 
sociated with magnetism and superconductivity and is a The Peierls instability towards a spatially broken sym- 

promising toy model for the electronic degrees of freedom metry state in one dimensional systems is caused by the 

of high-temperature superconductors. The most interest- competition between the energy of lattice distortions and 

ing properties occur mostly in underdoped samples, with the formation of gap at the Fermi level in the electronic 

electron densities close to half filling, where the system is spectrum. This instability can occur in two dimensions if 
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the structure of the Fermi surface has a strong nesting by 
a single vector. It is the case of the two-dimensional square 
lattice in the tight-binding approximation at half filling. 

Elastic Umklapp scattering with momentum transfer 
of (tt, 7r)[lJ across the Fermi Surface is allowed when it ex- 
tends to the Brillouin zone boundary at half filling. In this 
case, a small electron-lattice coupling will induce a lattice 
dimerization which is related with a periodic modulation 
of the bond hopping called in the literature, bond-order 
wave (BOW). This motivates the study of the typical pat- 
terns for elastic deformations with modes corresponding 
to the nesting vector (tTjTt). 

Although a finite frequency for the phonons would be 
important to study their influence in the mechanism of 
high Tc superconductivity, considering them as adiabatic 
could shed light on their influence on the different inho- 
mogeneous phases that have been observed, e.g. in the 
underdoped region. Besides, their role is expected to be 
important for the undoped parent compounds. 

Electron phonon coupling can lead to charge inhomo- 
geneities, such as stripes, as has been studied in [2\ within 
a spin-fermion model. However, this happens far from half- 
filling and for sufficiently strong diagonal coupling. 

On the other hand, the lattice distortions change the 
Fermi Surface (FS) shape and the FS deformation due to 
the presence of interactions is a central question within 
the breakdown of the Fermi liquid theory. [3l4l5l6j 

In the strong coupling regime one can find not only a 
deformed FS, but may even find a different topology. The 
FS of the non interacting Hubbard model is closed around 



in the Two-Dimensional Lattice-Hubbard model. 

the origin in the reciprocal space while in the interacting 
case it can be a surface closed around the point (tt, tt). It 
is in this regime where the role of the elastic distortions 
can take an important place in the high temperature su- 
perconductor phases [7] 

In this paper we study the stability of lattice distor- 
tions in the presence of Coulomb interactions U. The crit- 
ical values Uc where the distortions are suppressed are 
studied and a schematic phase diagram is drawn. The 
spectral properties are analized in the region of param- 
eter space where the elastic deformations are favored and 
the FS shape is constructed from the renormalized disper- 
sion. The FS shape was studied for the Hubbard model in 
the last years |8I9I10I11I12I13| , but the influence of lattice 
distortions has not been taken into account. The Polaron 
formation in the Holstein-Hubbard model was investigated 
in [14] by means the slave-boson saddle-point approxima- 
tion and the t — J model with electron-phonon interaction 
was studied recently [15] for small systems. We show that 
these distortions have an influence in the FS shape leading 
to an anisotropic FS in the S phase. 

2 Second order Perturbation theory. 
2.1 The model 

We examine the Hubbard model, coupled with a classical 
phonon field , which describes spin-i fermions in a two 
dimensional square lattice with nearest-neighbor interac- 
tions, in the presence of lattice distortions in the adiabatic 
approximation. This model represent a useful toy model to 
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describe the physic in the presence of vibrational modes . the lattice and U is the on site Hubbard interaction. It is 
The dependence of the hopping amplitudes tij is assumed convenient to work with the dimensionless parameters 



to be linear in the lattice distortions. 



























(T i=l j=l 












hH.c) (1) 


K 



where is the transfer integral between nearest neighbor 




(c) 



Fig. 1. a) Lattice shape for the modulations corresponding 
to pattern (a), b) Fermi surfaces for various densities p in the 
two-dimensional free fermion case. For densities p < 1 the FS 
is centered in the .T-point and for p > 1 it is centered in the 
Z-point. c) Diagrams which contribute to the self energy up 
to second order in U . The wavy lines represent the on site 
repulsive interaction and the solid lines the free propagators. 



sites in the absence of distortions, 7T.i.j,o- 



the number operator for electrons of spin a at site i, a is 
the electron-phonon coupling, K is the elastic constant of 



(2) 
(3) 
(4) 

and in the following we fix the energy scale setting = \. 

We first review briefly the free case {U = 0) and 5^'"^ = 
0, where: 



eo(k) — — 2(cosfca; +cosfcy). 



(5) 



The FS for different densities p is shown in Figure [T]-b, for 
half- filling {p — 1) being a diamond centered at the origin 
with vertices at ±(0, tt) and ±(7r,0). 

Let us consider the two possible alternation patterns 
consistent with the nesting vector Q = (tt, tt): 

SI, = i-iy+'s^o or si^ = i-iy+^s^o 

si, = o si = {-iy+^y, (6) 

namely patterns (a) and (b) respectively [16| . The two 
dimensional Peierls instability in the pure-hopping case 
was studied by Tang and Hirsch 17J founding that at half 
filling pattern (a) is favored. Using a Mean Field (MF) ap- 
proach we can see that for a wide range of values of U this 
pattern survives up to a critical value Uc- First we review 
the MF results and later, by using second order perturba- 
tion theory we reexamine the stability of the distortions 
in the present model. 

When one considers a non trivial Hubbard local inter- 
action, it is easy to see, using Hartree-Fock [TB], that for 
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U < Uc distortions (a) are favored at half filling and for 
U > Uc the distortions are quickly suppressed. 

Treating the local Coulomb repulsion in the MF ap- 
proximation, 

- K,j,T)Kj,i)) (7) 




0.05 0.1 0.15 



0.25 0.3 



the hamiltonian becomes quadratic and it is straightfor- 
wardly diagonalized in reciprocal space. 

The expectation value {ui^j) is assumed uniform, and 
the electron density with a given spin can be assumed 
as (riij^a) = I + ^{ — ly^-'m pj], m being the staggered 
magnetization. Within this approximations the electronic 
spectra for the two patterns (a) and (b) are 



+ 4(OOB + COB ky)'' + 4S 



2 4- sin2 I 



[y2^2 

e('''(k)± = ±2t, h 4(ooB fc^ + COB fca)2 + 4a2(Bi„ fc^ + ai„fcj,)2 



The ground state energies for the two patterns are 



with r = 1 for pattern (a) and r = 2 for pattern (b). 

Minimizing the energy with respect to S and m for 
various values of U and A we can easily see that pattern 
(a) has lower energy than pattern (b) and for a huge range 
of U the values of S^jp for minimum energy are almost 
constant. 

For example for A — 0.5 we have for pattern (a) ~ 
0.24 and for pattern (b) Sl^p ~ 0.13. The value of the 
interaction where the distortions are suppressed is Uc 2. 

In the presence of distortions that follow pattern (a) 
the lattice changes as we show in Figure [T]-a. The unit 
cell is doubled and the Brillouin Zone is reduced to half. 
It is straightforward to find the electronic spectra for the 



Fig. 2. Energy calculated for both patterns (a) stars and (b) 
squares as a function of S for U — 0.4, A = 1/2 and m = 0. 
Clearly pattern (a) has lower energy. On the other hand it is 
easy to see that the value of 5 in the minimum energy is almost 
constant in the range of U between and ~ 2. In the Inset the 
difference of energy between the two patterns is shown as a 
function of U, this difference falls off to zero at (7c ~ 2 where 
the distortion 5 its suppressed. 

deformed free case {U = 0). Working in the first Brillouin 
zone (BZ) we have two single particle bands 



£o(k)± — ±2tW (cos kx + cos ky)^ + sin k^ 



for pattern (a) and 



£o(k)± = i'2t^ (cos kx + cos ky)'^ + (5^(sin'^ kx + sin^ ky). 

for pattern (b). In the following we will work in the ex- 
tended Brillouin zone formed by a square centered in the 
point r = (0,0) with vertices in ±Z and ±Z where Z — 
(tt, tt) and Z = (vr, — vr). The unperturbed dispersion in 
the extended BZ is 



;(k) = 2(1-2 ei(j)(k)))^jj(k)2 + 5^ sin^ 



(10) 



£(k) = 2(1-2 e(r,(k)))^jj(k)2 + 52(sin2 + sin^ ky) (11) 

for patterns (a) and (b) respectively, where we have used 
?7(k) ~ cos kx + cos ky and 9 is the step function. In order 
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to fix the notation we present here other important points 
used in this work hke X — {Tr,0), Y — (0,7r) and M — 

( IL 1L\ 

V 2 ' 2 

2.2 Second order self energy 

The contributions to the self energy are given by the dia- 
grams shown in Figure. [T]-c, these contributions were cal- 
culated using ordinary perturbation theory and were cal- 
culated slightly away half filling where divergences can 
appear. To deviate sUghtly half- filling are not significant 
changes in patterns of deformations in the mean field treat- 
ments and previous work showing that there is no appear- 
ance of stripes phase for fillings close to one [5]. 

Since the partners in the scattering processes have op- 
posite spins, other possible diagrams with two Coulomb 
lines are absent. 

The contributions of these diagrams are given by 

-OO JbZ 



2 



(12) 



where p is the electronic density of the system and 
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Fig. 3. Self Energy in M, F ,X and Z points at p = 0.995, 
The solid line correspond to the real part and the dotted one 
to the imaginary part of the Self Energy. 

and does not contribute to the FS deformation because 
it can be absorbed by a shift in the chemical potential 
Sjii — Un/2 to keep the density fixed. The second order 
contribution to the self energy 



L'^ k"^" ^ ^ ^ Ck-k' - Ck'-k" + ^k" 
Ck-k')^(~Ck'-k" )^(Ck" 



(Stt)^ J _aa J-oo J BZ J BZ 

X g(°) (k', (^')g'"' (k", ^'Og'"' (k + k' - k", + 1^' - 



where 



G(")(k,c^) = 



1 



(14) 



W - Ck + «7k 

with = f o(k) — /i, 7k = 7 sign(^k) and [i is the chemical 
potential. 

To first order, the contribution of Z'^^) is a real con- 
stant (k-independent) that shifts the dispersion relation 



— 4k-k' — c;k'-k" + 4k" / 

is k-dependent producing a renormalized dispersion and 
leads to a FS deformation. This term is computed nu- 
itli§^ically in what follows using a Fast Fourier Transform 
L<('f)F,T) for finite size systems. The momenta in the Bril- 
louin zone are discrete and defined by kx,y — —Tr+Ak^n^^y- 
1) and Ak = ^, n,,y = 1, 2, L 

The retarded self energy in space-time representation 



IS 



sift (x, t) = TuT^ (k, iO+ 

_ 1 ^-^ duj 

~ 7.^ 



(16) 

e*('^-^-'^*)Z'(2)(k,w + iO+) 
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Fig. 4. Renormalized dispersion relation corresponding to 
5 = 0.24 and (7 = 1.8 in the T ^ X ^ ^ T line. The 
squares correspond to the renormalized dispersion calculated 
by means of Eq. (|21|l and the solid line corresponds to the 
free case. We can see that the Coulomb interaction reduces the 
bandwidth. In the ^-point some high energy solutions appear. 

where the Fourier transform is defined by 

^k[5(k)](x) = ^^e^''-g(k) (17) 

k 

/OO 7 
^e-^-'giio) (18) 

and the inverse transformation is denoted by !F^^ and 
respectively 

The frequency Fourier transform is straightforward 



^ z(k - k\ t)z(k' - k",t)iu*(k",t) 

-J,k" 

+ tu(k-k',t)tu(k'-k",t)z*(k",t)] (19) 



where we have used the short-hand notation 
z(k,t) = e-^«''*0(-a) 

The transformation in k is the Fourier transform of a con- 
volution and can be written as product of the Fourier 



transforms 

rfj (x, t) = -ie{t)— [i" (x, t)w (x, t) + (x, t)s* (x, t)] 

with 2(x, t) = jFk[2:(k, t)] and w(x, t) = J-k[i(;(k, t)] 

Clearly i(x,i) = (cx(t)4(0)} and^I;(x,i) = (4(0)cx(0) 

and the causal free propagator in space-time is given by 

G°(x, t) = -i(6l(<)5(x, t) - 0{~t)w{yi, t)) 

The functions z(x, t) and w(x, t) can be calculated for 

finite size systems with a standard FFT algorithm. Then, 

(2) 

replacing the result in the equation for S^gt (x, t) , we can 
obtain the self energy in momentum space by means of 
the inverse Fourier transform 



The real and imaginary parts of the self energy are 
shown in Figure[3]for M, F, X and Z points at p = 0.995. 
The shape of the Self Energy is very similar to the one 
found for the Hubbard model [19] but some features are 
(5-dependent. For example, in Figure[3]the real part of U in 
the point X has a clear linear behavior in a wide interval 
around lu = 0. The range of to where this occurs is larger 
for higher values of S. The symmetries of the Hubbard 
model are preserved. We can see in Figure [3] that if one 
changes lo —u in the plot for the imaginary part of S 
at point r, the one corresponding to point Z is recovered. 
In the plots for the points F and Z it can be seen that 
?fte{S{k,uj)} = -5Rc{Z'(k + Q,-w)}, with Q = (tt.tt). 

The interacting Green function [5D] can be calculated 
up to second order 



G(k,cj + iO+) 



w + iO+ - Ck - ^(k, LU + iO+) 



(20) 
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and the low energy excitations can be determined from 
the equation 



(21) 



For each point in the Brillouin zone, the last equa- 
tion gives the renormalized dispersion as a function of the 
chemical potential 



(22) 



In Figure [5] we show the renormalized dispersion for 
U — 1.8 in the path shown in Figure [TJ We can see that 
the Coulomb interaction reduces the bandwidth. 

For values of the momentum away from the FS, in 
particular for high values of U , Eq. (j2ip can give more 
solutions corresponding to higher energy excitations. This 
is particularly clear in the points 7^ and Z where two new 
bands appear for U > 1. This is reminiscent of the Hub- 
bard bands. Figure 3] shows the high energy solutions in 
the Z-point. We chose to show a result for a high value of 
U because the effects of renormalization are more evident 
and then the high energy solutions arc visible. 

2.3 Stability of lattice distortions. 

The Ground State energy per site was calculated for sev- 
eral values of U and A 
2 



E = 



J d'kS^{S)9{-S^{S))+r^ (23) 



where r — 1 for pattern (a) and r = 2 for pattern (6) 

Keeping fixed the filling, we can calculate the energy 
for various values of S and A as a function of U. In this way 
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Fig. 5. Energy per site vs U for several values of 5 correspond- 
ing to A = 0.5 and p — 0.995. For Uc ~ 1.21 there is a crossing 
between the curves and for U > Uc the minimum energy cor- 
responds to 5 = 0. In the inset we show the optimal value 5* 
that minimizes the energy as a function of U. For U > Uc the 
distortions are suppressed and for < U < Uc the parameter 5 
takes an optimal value S* for most of the values of U and there 
is one narrow region near Uc where S goes quickly to zero. 

we can study the stability of the lattice distortions when 
increasing the Coulomb interaction U, beyond the Mean 
Field approximation. In the following we restricted our- 
selves to study pattern (a) . The values of S that minimize 
the total energy depend of A and when A ^ 0, the distor- 
tions are suppressed. In our calculations the deformations 
are suppressed for values below some A finite since the 
calculations were made for finite size systems (in general 
for systems of 120 x 120 sites, but the spectral properties 
are not severely changed for bigger sizes. ) In Figure [5] we 
show the energy as a function of [/ for A = 0.5 and sev- 
eral values of 6. We can see that for U < Uc the curves 
have a minimum for a non trivial value of S and that for 



u 
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spectral properties, we labeled this phase simply as 5 = 
or non-elastic phase. In the region below the critical curve, 
where the elastic deformations becomes stable, the phase 
is labeled by (5 7^ or clastic phase. For points close and 
below the critical curve there is a small region where the 
5- values rise from to (5*, as can be seen in the inset of 
Figure [S] 

For fixed values of the Coulomb interaction the electron- 
lattice constant A governs the phase of the system, and for 
this reason it could be interesting to compare the prop- 
erties of the system at fixed U ioi S — and 6 = 5*. 
In the following we show some of the spectral properties 
calculated for the model in the elastic phase. 





u = u. 


5 = 









0.2 



0.4 0.6 . 

A 



Fig. 6. Schematic phase diagram. The curve corresponds to 
U = U,. 



U > Uc the minimum energy is reached for 6 = 0. The 
5*-value for minimum energy is plotted as a function of U 
in the Inset of Figure [5] where we can see clearly that S* 
stays almost constant for values of U smaller than Uc and 
it tends quickly to zero in a narrow region U ^ Uc. The 
critical value Uc that we found in our approach is lower 
than the one found in [15] by means of a MF approach, 
but the corresponding values of 6* are very similar. 

In Figure [6] we show schematically the regions in the 
(A, U) space where 6 takes non trivial values, this region is 
labelled as (5 ^ 0. The curve in Figure [S] shows the values 
found for Uc- For each value of A we have 6 = for the 
[/-values above the critical curve. The two dimensional 
Hubbard model with [/ > at half filling has AFM long 
range order, however, to depart slightly from this point we 
can not ensure that this order is maintained in the sys- 
tem without properly studying the AFM order parameter, 
but as the main subject of this work is the study of the 
stability of elastic deformations and their influence in the 



2.4 Spectral properties in the presence of elastic 
distortions . 

The single-particle spectral function ylk(ti') = — iSmG(k, a;) 
is calculated in the second order approximation. In Figure 
[7] we show the spectral function along the lines F ^ Z, 
Y ^ Z, X ^ Z and r ^ X for U = 1 in the 5 ^ 
phase. The behavior of the spectral function is seen to 
be different from the one found in the region with 6 = 
|19j . The peaks in the spectral function for this case are 
broader and asymmetric for k far from the FS, while for 
the case with distortions the symmetry is preserved and 
the peaks are narrower. In the 6 = case at higher ener- 
gies the quasiparticle peaks decay very slowly while in the 
6 case the peaks decay more quickly. The splitting 
in the spectral function seen for (5 = is absent in the 
(5^0 case. In the presence of distortions we can see the 
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Fig. 7. Spectral function for the lines F ^ Z,Y ^ Z, X ^ Z 
and r —* X. The line in Fig. (a) is crossing the M point where 
there is no spectral peak as a result of the presence of a gap in 

Fig. 8. a) Density of states p(uo) in the case with 5 = for 

the quasiparticle dispersion. 

U = 0.8, f/ = 1, (7 = 1.2 and U = 1.3. The interaction gives a 
absence of a spectral peak for some points around w = transfer of spectral weight from low to high energies, b) Density 
which shows the existence of a gap at these points in the of states p{uj) in the presence of elastic distortions for U = 0.8, 
renormalized dispersion. There are no possible low-energy U = 1, U = 1.2 and U = 1.3. 



excitations with r; « at these points of the BZ. 

The density of states was calculated for fixed values of 
the chemical potential and the Coulomb interaction U . 



give rise to the typical two-peaks situation corresponding 
to the infinite U limit. 



(24) 



In Figure [8] the renormalized density of states is shown for 
both phases at p = 0.995 with U = 0.8, U = 1, U = 1.2 
and U — 1.3. The interaction gives a transfer of spectral 
weight from low to high energies in both cases. 

In the absence of distortions the weight of the logarith- 
mic singularity which characterizes the free system is re- 
duced. Similar characteristics of g{e) are obtained for the 
infinite-dimensional Hubbard model [21j . Large U values 



2.5 Fermi Surface in the presence of distortions 

The solutions of w — — 5ReZ'(k, w) = determine the 
renormalized dispersion S'k, the points where S'^ — de- 
fine the interacting FS. When we analyze the Hubbard 
Model without elastic distortions a point with S'k < is 
inside of the Fermi Area and if ^E'k > it is outside. In the 
present case we need to be careful with this interpretation; 
a point with S'k < can be inside the Fermi Area or at 
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3 Discussion and Summary 



ky 





14 — 
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Fig. 9. Interacting Fermi Surface for U — 0.8, fi = —0.1 and 
5 = 0.24. The red circles mark the zones where the surface 
stretches in the fc^i-direction and the blue arrows the zones 
where the surface is contracted in the y-direction. 



the FS depending on whether the point coincides with the 
gap in the renormahzed dispersion relation. 

The interacting FS for S = 0.24 and /i = —0.1 is shown 
in Figure m The symbols correspond to U = 0.8 and the 
solid line to U — 0. We can see that the effect of the in- 
teraction leads to an anisotropic surface resembling a ne- 
matic phase FS. The points in the interacting FS evolve 
so that the point on the fc^j-axis comes closer to the X 
point and the point on the fcj,-axis moves away from the 
Y point when the interaction is increased. We do not see 
any change in FS topology. This result is consistent with 
earlier works 12J and shows that interactions do not mod- 
ify the FS topology within the perturbatively controlled 
weak coupling regime. 



The Hubbard model on the square lattice in the weak 
coupling regime in the presence of lattice distortions that 
follow Peierls-like patterns was investigated. 

Using second order perturbation theory, several spec- 
tral properties are calculated and compared with the Hub- 
bard model in the absence of lattice distortions. The sta- 
bility of the distortions was analyzed as a function of U 
and A finding that the Coulomb interaction suppresses the 
lattice distortions for values of Uc smaller than previous 
MF results and a schematic phase diagram is presented. 

The results show that the Interacting Fermi Surface 
is anisotropic in the presence of distortions, even in the 
weak coupling regime. The Fermi Surface topology does 
not change in any of the two phases. This result is con- 
sistent and complementary with earlier results [12] where 
the stability of the FS topology in the absence of distor- 
tion has been analyzed previously by Metzner et al in the 
weak coupling regime . The results presented in this paper 
show a similar behavior in the presence of distortions. 
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